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Abstract 

The standard Born Oppenheimer theory does not give an accurate 
description of the wave function near points of level crossing. We give 
such a description near an isotropic conic crossing, for energies close to 
the crossing energy. This leads to the study of two coupled second or- 
der ordinary differential equations whose solution is described in terms 
of the generalized hyper geometric functions of the kind 0-^3(1 cl^ b, c; z). 
We find that, at low angular momenta, the mixing due to crossing is 
surprisingly large, scaling like fJ^^^^, where fi is the electron to nuclear 
mass ratio. 

1 Introduction 

In 1927, in a landmark paper. Born and Oppenheimer P| paved the way to 
applying quantum mechanics to molecular spectra. In their paper they in- 
troduced an approximation that greatly simplified the treatment of quantum 
mechanical spectral problems in which the particles can be divided into heavy 
and light. Molecules are an example since the nuclei are much heavier than 
the electrons. We shall denote by ^ the small parameter of the theory. In 
molecules /i ~ 10~^, the electron to nucleon mass ratio. Since the light par- 
ticles are associated with the fast degrees of freedom and the heavy particles 
with the slow degrees of freedom, the Born Oppenheimer approximation is 
related to the adiabatic approximation |^, |l^]. At the same time, the Born 
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Oppenheimer method can be viewed as a generalized semiclassical approxi- 
mation where the small parameter plays the role of h^. This is a reflection 
of the fact that the electrons and nuclei also live on different spatial scales, 
the electronic wave function is spread and is far from the semiclassical limit, 
while the nuclear wave function is tight and close to semiclassical. 

The procedure put forward by Born and Oppenheimer is to flrst solve 
the electronic spectral problem with flxed nuclei, and view the nuclear co- 
ordinates as parameters. To the leading order in /z, and far from crossings, 
the heavy degrees of freedom in Born Oppenheimer theory are described by 
a (scalar) Schrodinger operator in the semiclassical limit with = n and 
where the electronic energy surface serves as a potential. 

Born Oppenheimer theory developed into two distinct directions. The 
main direction has been the application to various systems and the devel- 
opment of effective and accurate methods of calculations ^, [TD|, |TT]]. The 
second direction has been the development of the theory as a tool of rigorous 
spectral theory |^ H, || . Our work falls into the flrst class. 

Points where electronic energy surfaces cross are singular points of the 
Born Oppenheimer theory. In certain cases, these points can affect spectral 
properties P]. In this work we focus on the behaviour of eigenfunctions near 
crossing. There is surprising little that is known about this. It is not known 
if the function has flnite values at the crossing; how the amplitude of the 
wave function near the crossing scales with fi. Since crossing controls the 
mixing of electronic levels, the knowledge of the wave function near crossing 
is important. Of course, these questions are interesting in the case that the 
crossing lies in the classically allowed region. 

In this work we address these issues for conical (i.e. linear) crossing of 
two levels where the Born Oppenheimer problem reduces to two coupled 
Schrodinger equations. In the isotropic case the analysis reduces further to 
the study of two coupled, second order, ordinary differential equations. We 
obtain the nuclear wave function analytically, to leading order in fi, close 
to the crossing. It is related to the generalized hypergeometric functions of 
the kind 0-P3. This function takes the ordinary Born Oppenheimer nuclear 
wave functions, which are a good approximation far from the crossing, all the 
way to the crossing, where the nuclear wave function mixes the two electronic 
levels. We flnd that the nuclear wave function with total angular momentum^ 

^We do not actually mean here the physical angular momentum, but a quantum number 
reminiscent to it, see Eq. (|l2|). 
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m = ±2 is nonzero at the crossing point. Moreover, for low momenta, the 
wave function has a large amplitude near the crossing, of order yu~^/^. We 
find appreciable mixing of the two levels at distances that are smaller than 
and the total weight that is mixed between levels scales like ij}/^ which 
is remarkably large. 

2 The Born Oppenheimer approximation 

This section is a brief introduction to the basic and elementary elements of 
Born Oppenheimer theory. 

2.1 The basic model 

A prototype of the Born Oppenheimer problem, and the one we study here 
is§: 

H = -fiA, + Heix) (1) 

where /i is a small parameter. He{x) is an operator valued function that acts 
in the Hilbert space of the light degrees of freedom. Since the eigenvalues 
of He{x) repel |T5[, one does not expect crossing in the case of one heavy 
coordinate x. If He{x) is time reversal invariant, then stable crossing will 
occur if there are two heavy coordinates. We therefore assume that a; G M^. 
For the sake of simplicity we have taken identical masses for the two heavy 
degrees of freedom^. Two degrees of freedom is the simplest case that is still 
rich enough to cover the phenomena we are interested in. 

There are several ways to motivate H. The most direct is to think of H 
as a phenomenological quantization of the molecular vibration. For example, 
in the case of molecular trimers the two heavy modes are the antisymmetric 
stretching and bending of the molecule, see fig |^. 

Alternatively, one can start with the Schrodinger equation for a molecule, 
which indeed has the form of Eq. (|T|) where Hf.{x) includes the Coulumb po- 
tential of the nuclei and the electrons and the electronic kinetic energy. Often, 
and this is the case in molecules, if, is invariant under Gallilean translation 
and rigid rotation. The symmetry gives three quantum numbers {P, J, m} 
and the spectral analysis of -ffp,j,m. is now restricted to the "internal" nuclear 
coordinates. In |p one can find a detailed description of this procedure for a 

^There is no loss here, for by scaling some of the x-directions this can always be achieved 
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Figure 1: Two degenerate modes in triangular molecules. 



triatomic molecule. Hpj^m has a more complicated expression than H: Fix- 
ing the center of mass at the origin and restricting to an angular momentum 
subspace replaces the kinetic energy — /i A by a more general quadratic func- 
tion of the of the momenta. However, locally near the crossing this expression 
reduces to § Eq. (§. 

We shall assume that He{x) has discrete spectrum and has smooth depen- 
dence on the coordinates^. In addition, since we shall use the Born Oppen- 
heimer theory as a calculational tool, rather than a tool of spectral analysis, 
we shall assume that the problem has benign qualitative spectral features. 
For example, we shall assume that the spectrum of H in the energy range of 
interest is discrete, and that the associated wave functions are localized in 
space in the classically allowed region, and that this region is connected. Sub- 
tleties associated with tunneling and other exponentially small phenomena 
will not concern us here. 

Since the small parameter /x multiplies the leading derivative in the x 
variable, the Born Oppenheimer problem is a version of semiclassical problem 
where the operator Hf.{x) replaces the scalar potential V{x). 

In molecules there is a second small parameter, 1/c ~ 1/137, which gov- 
erns relativistic effects. In particular, spin-orbit interactions, are of order 
1/c^. The lowest order of Born Oppenheimer theory gives an energy scale of 
^//I which is comparable to 1/c but is much larger than 1/c^. It is therefore 
consistent when discussing Born Oppenheimer to leading order to disregard 
spin-orbit. This is also the reason why we shall not go beyond the leading 

•^Both assumptions are not realistic; The electronic Hamiltonian has continuous spec- 
trum at high energy, and because of Coulombic singularities there is no smoothness in 
the X dependence. Fortunately, both problems are, by now, well understood and may be 
viewed as a technical complication that, for our purposes, can be left out. 
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order. 



2.2 Partial Diagonalization 

The starting point of the Born Oppenheimer method is to consider H in the 
basis that diagonahzes the fast (electronic) degrees of freedom. 

We assume that the electronic Hamiltonian is real, which is the case in the 
absence of external magnetic fields. Let 0{x) be the orthogonal transforma- 
tion that diagonalizes He{x). If He{x) has simple (non-degenerate) spectrum 
in the vicinity of x then 0{x) is uniquely determined up to multiplication 
by a diagonal matrix with ±1 on the diagonal. Locally, one can choose 0{x) 
so that it inherits the smoothness properties of Hf,{x) |T^. It follows that in 
the basis that diagonalizes He{x), Eq. (|I|) takes the form 

O^HO = fi{-iV^ - A{x)f + E{x), (2) 

where E{x) is a diagonal matrix whose entries Ej{x) are the electronic energy 
surfaces and A{x) = iO''{x)'VO{x) is a (matrix) gauge field. Since 0{x) is 
real, A is antisymmetric and hermitian. A is responsible to the coupling 
between electronic energy levels. 

One can associate to a crossing point n indices, {ii, . . . ,in}, where n = 
dim He and ij = ±1. Let us take a closed curve around a crossing point. 
After such a cycle 0{x) must return to itself up to up to multiplication by 
a diagonal matrix with entries ±1 on the diagonal. These entries are the 
indices of the crossing. It is known that for a conic crossing between the 
j — th and j + 1 eigenvalues of H^^x) give ij = ij^i = — 1 and all other indices 
are, of course, 

For conic (linear) crossing 0{x) flips signs on a circle of radius 0. 
Therefore its gradient must be of order l/|x|. This makes A of order l/\x\. 
It follows that the coupling between electronic states diverges like a simple 
pole near crossing. 

This completes the local description of the theory. The problem also 
has an interesting global aspect. Since He{x) is a real symmetric matrix, 
the Wigner von Neumann crossing rule ||15[| says that iJe( G M has, 

generically, isolated crossing points. As a consequence, with points of crossing 
being removed the plane becomes multiply connected (see fig. |) 

We can now describe the boundary conditions associated with Eq. (^. 
The general case of several crossing points can be complicated but in the 



5 



Figure 2: With points of crossing removed the plane becomes multiply connected. 



case of at most one point of crossing the situation is simple. In that case, 
cut the plane from the crossing point to infinity. On the cut plane 0{x) is 
uniquely defined in a continuous way. Then the boundary condition on the 
j-th component of the wave function associated with Eq. is periodic or 
anti-periodic according to the index ij. 

2.2.1 Born Oppenheimer theory near a non-degenerate minimum 

Let x = be a minimizer of an electronic energy surface. Pick the origin so 
that the minimum is at zero energy (see Fig. Upon scaling, x = n^^^^, 



E{x) 



X 

Figure 3: A potential energy surface near a non-degenerate minimum, 
the Born Oppenheimer operator assumes the form 

(3) 



((-^Vg - e A{ei)f + E{ei) 



/i 



The (scaled, matrix) potential energy is 0(1) for the electronic energy surface 
near the minimum, and has gaps of order to (scaled) "excited electronic 
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states". Suppose first that there is no crossing. Then the couphng between 
electronic levels is small, e A{e^) = 0{ij}^^), and a perturbation argument 
shows that the effect on eigenvalues is of order 0(/i) (in unsealed energy) and 
of order 0(yU^/^) for the wave function. It follows that the spectral analysis 
in an energy interval of order o(l) near the minimum, reduces to an ordinary 
Schrodinger equation (not matrix valued) with no vector potential (since A 
vanishes on the diagonal). This accounts for o(l//i) eigenvalues near the 
minimum in two dimensions. 

Now, if there is a point of crossing, there are two possibilities. The first, 
and simplest, is that the crossing lies in the classically forbidden region, 
see Fig ^. Then the divergent vector potential is harmless, since the wave 




Figure 4: For a non-degenerate minimum, when the crossing falls in a classically 
forbidden region, the cut to infinity can be pushed to the classically forbidden 
zone, and one can forget about the crosing to leading order. 

function is exp(— -^) near the crossing. The cut to infinity can be pushed 
to the classically forbidden zone, so the difference between periodic or anti- 
periodic boundary conditions is exponentially small, and one can forget about 
the crossing altogether. If the crossing point lies in the classically allowed 
region, it couples two nuclear Scrodinger equations and ruins the traditional 
Born Oppenheimer approximation. 

2.2.2 Born Oppenheimer theory near a degenerate minimizer 

It can happen that the crossing must be taken into account although it lies 
at the classically forbidden zone. This happens when the cut can not be 
pushed to the classically forbidden zone. This is the case, for example, when 
the curve 7 G is a minimizer of an electronic energy surface with zero 
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energy Q on 7. If 7 encircles a crossing (see Fig. ^ then the cut necessarily 




Figure 5: A degenerate minimum around a crossing. The solid circle represents 
7. Although the crossing falls in the classically forbidden zone, it should be taken 
into account. 

intersects 7 and since 7 lies in the classically accessible region, the wave 
function is large there and it matters if one imposes periodic or anti-periodic 
boundary conditions. 



3 Born Oppenheimer theory near crossing 

Consider now the spectral problem near crossing energy of two electronic 
energy surfaces. Let us set the crossing energy at and assume that the 
crossing is conic. This is the generic situation. 

Upon scaling, x = e^, e = /i^''^, the Born Oppenheimer operator assumes 
the form 

e(^{-tV^-eA{eOY+lE{eOy 6 = //3. (4) 

The scaling increases the gaps between the electronic energy surfaces and 
decreases the coupling of the crossing pair to other levels since eAij{e^) = 
0{ij}/^) when i belongs to the pair and j to other levels. On the other hand, 
the crossing pair remains coupled, because the Coulombic singularity near 
the crossing says that ey4j^j/(e^) = O (^j^jj? which is large when ^ is small. 
We see that a spectral problem in an interval of order o(l) near the crossing 
^Similar arguments apply if the energy associated with 7 is sufficiently small 
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energy, reduces to a problem where He (x) is a 2 x 2 matrix up to an error of 
order /x in the eigenvalues and of order jj?/^ in the eigenf unctions. 

Our aim, in this work, is to describe the wave functions for states, located 
at an interval of width that is much smaller than /x^/^ near the crossing. 
Even though this is a small interval, it has lots of eigenvalues: By Weyl's 
rule there are many states, of order /x"^''^, in an interval of width /i^^^, in two 
dimensions. One can expect to find many states in the interval in question. 

Close to the crossing He{x) can be expanded in terms of the Pauli matrices 
and the unit matrix. We assume that asymptotically close to the crossing 
point He{x) is isotropic and conicQ: 



He{x) = a-x + 0{x^). (5) 



We choose (Ti,2 realQ 



^^=(? J)' ^^=(J -l)- 

0{x'^) accounts for the behaviour far from the crossing which is not universal, 
or isotropic. 

How restrictive is the assumption that near the crossing He{x) is isotropic? 
Molecules are never isotropic, although some may be approximately so. Nev- 
ertheless, isotropic conic intersections are, in fact, common. This is a known 
phenomenon in group theory: Discrete symmetry can force full continuous 
symmetry on tensors of finite rank [^. For example, as shown by @, the 



D3, symmetry of trimers forces isotropy of the conies. 

It follows from the above that if one is interested in the local behavior 
of eigenfunctions for eigenvalues that lie in an energy range that is small 
compared to fj.^^^ and in a spatial neighborhood of the crossing, |x| << 1, 
the eigenfunctions satisfy, to leading order in /i, a canonical system of partial 
differential equations, 

(-/iA,, + a-a;)^ = /ii/3(- A^^ + a-i)^ = Q. (7) 

The rotational symmetry allows us to reduce the system Eq. (|^), to a system 
of linear, ordinary differential equation parameterized by angular momentum 



""Reality, isotropy and linearity would allow for an additional overall scale factor in Hf. 
This scale can be absorbed in a redefinition of /i. 
^This is unconventional, but convenient. 
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m G Z + I (see section ^?T| ): 



For fixed m the space of solutions of the ODE is four dimensional. We shall 
see that there is a one dimensional subspace of solutions that is well behaved 
near the origin p = 0, and near infinity. As we shall explain, functions in 
this space describe the asymptotic behavior of eigenfunctions with energies 
near the crossing and spatially close to it. The explicit expression for these 
solutions, is described below. 

3.1 The main result 

We now describe our main result: 

Theorem 3.1 For m G Z + and m > ^, let 

f I „\ _ ^ " l_ ( O-' 3li 3) 2 ^ 3 ' 6 ^ 3 ' e'*^' I 

^ I3M I2 + U + V e^0^3l, 3' 2 + T' 6 + T' F V 

r(i)r(| + f)r(| + f) V oF,{-M + %l + t4) ) 

then 

• ^m{p), is a solution of the system of differential equations, Eq. (Qj. 

• For small p 

(0)^ i r(i)r(i+f )6i ^ /g>) 

l(2 + y)b 5 Y r(i)r(|+f)/ 

In particular, J^i(O) 7^ 0. 

• For large p 

- pj)^ (fp"^ - ^ (f + 1) ) ( -1 ) 

10 
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Figure 6: The components of J-m{p) for m = ^ and their asymptotic form given by 
([l0|). The function is compared with the WKB Born Oppenheimer wave function. 

Theorem 3.2 For m E Z + ^, and m>\, let 

f pie 12 Ap~ie/2 \ 
^^{r,e) = i^-'"e-' i^^,,/, j T^{yr'l\) (11) 

then 

• is a solution of the system of PDF, Eq. (Qj where x = {r cos 6, r sin 6) . 

• The m-th component of an eigenfunction of Eq. near crossing, i.e. 
with eigenvalue << fj,^^^ and for x = (r cos 6*, r sin 6*), r « 1 is, to 
leading order, proportional to \l/m('"7^)- 

• The amplitude of'^m is independent of jj, in the region r» jjL^/^ . 

• Near the crossing, r ^ jJi^^^ , the amplitude of the wave function = 
0(/i-V4) 

Remark 3.1 The most interesting aspect of the solution is that the wave 
function has large amplitude at the crossing region in the limit of small /x. 
As we shall see, this result follows from arguments that do not rely on the 
explicit form of the solution, hut do depend on the fact that in the non-mixing 
region, r» fj.^^^ the solution has a WKB form in the radial direction. 



11 



Remcirk 3.2 The function does not describe the behavior of wave func- 
tions in the far zone where r > 1. The behaviour in the far zone is not 
universal and depends on the details of the electronic energy surface . The 
far zone is described by standard Born Oppenheimer, so gives comple- 
mentary information. 

Remark 3.3 There are a (/i^^/'^) states in the relevant energy intervale [fi^^^) 
near the crossing disregarding m. For a given m there are only fi'^^^ levels in 
this interval. This says that near the crossing m is bounded by order fj,"^^^. 

Remark 3.4 We shall actually only prove the theorem in the special case 
that He{x) is rotationally symmetric. The symmetry decouples channels with 
different angular momenta. We believe that mixing of angular momenta in 
the far zone is only a technical complication and that the result also holds 
without rotational symmetry in the far zone. 



4 The Rotationally Symmetric Case 

In the foUowing we describe a derivation of the main result for a Born Op- 
penheimer model that is rotationally symmetric. With rotational symmetry 
we can reduce the spectral problem of a PDE to a spectral problem of an 
ODE. No real molecule is rotationally symmetric and the general case leads to 
mixing of m channels. We believe that this is only a technical complication. 

We require invariance of Hf.{x) under infinitesimal rotations in the nuclear 
and electronic Hilbert spaces. Such a rotation is generated by 

1 . d . d I 

Js-Ls + -as = -zx,— + ^x,— + -as, (12) 

with 

/ i 
""'^[-i 

The Ls part generates SO{2) rotations in the nuclear Hilbert space {xi — X2 
plane), whereas the |o"3 part generates a rotation in the electronic Hilbert 
space. Js does not have the meaning of total angular momentum since the 
Pauli matrices do not represent spin. Isotropy means that Js commutes with 
Heix): 

0^[Js,He{x)] 
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The most general form of Hf.{x) for a two level system that is rotationally 
symmetric and real is: 

He{x) = Qo{r) + Qi{r){x-a) + Q2{r){xxa), r = \x\ (13) 

An additional Q3{r)a^ term in (^) is allowed by rotational invariance but it 
is forbidden by time reversal symmetry, since Us is imaginary. 
The energy surfaces of He{x) are equal to 

Ei{r) = Qo{r) ± g(r); q{r) = r ^ Ql{r) + Ql{r) (14) 

We shall assume that (5o,i,2 are smooth functions of r and that (5i(0) = 1 
while Qo,2('") = O(r^) for small r. This gives conic intersection at zero with 
E±{r) = ±r, see Fig. ^. 



0.4 
0.2 ~ 




Figure 7: The electronic energy surfaces of He{x) 



4.1 The radial Hamiltonian 

The spectral subspace of J3, with eigenvalue m is spanned by 

m must be half odd integer for 61^2 to be univalued, i.e. m G Z + |. 
Since 

-Aei = ^^ei, -A 62 = ^^62, 



^(32(0) can be always set to zero by an appropriate rotation in the x space only, i. e. by 
an appropriate choice of the heavy coordinates. 
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x-aei = re2, x-ae2 = rei, 
X X a ei = —ir 62, x x cr 62 = ir Ci, (16) 

in terms of the basis {ei, 62} we obtain for the radial equation: 

H{m, ;x) = -/X + - + H,{r, m) (17) 

with 

He{r,fi,m) = Qo{r) + rQi{r)ai + rQ2{r)a2. - ^ {ma2 - m^) (18) 

Scahng r = yU^/^p we get Eq. (|^), to leading order in /i, for < p << p^^/^. 

Remark 4.1 The radial Hamiltonian, He{m, p, r), actually has no level cross- 
ing. This, by itself, does not ameliorate the mixing of the two levels for now 
the gap in the spectrum of He{m, fx,r) is of order jj}/^ . The smallness of this 
gap leads to mixing of the electronic levels. 
We shall restrict ourselves to m > 0. Since 

aiH*{r,fi,m)ai = He{r,ii,-m), (19) 

H{m,^) and H{—m,ii) are isospectral and the radial part of the function 
with — m can be obtained from the one with +m by interchanging upper and 
lower components and taking complex conjugates. 



4.2 The Indicial equation 

The origin, p = is regular-singular point pOf of the equation. Substituting 

(^;^)(l + 0(p)) (20) 

into (^ we obtain the roots a = ±(m — |) and /5 = ±(m + |). Equation (|[) 
therefore has four linearly independent solutions, which asymptotically near 
the origin, behave like 



p™-2 \ / p-™+2 \ . / \ / 



p™+2 ; ' \ p - 2 



(21) 
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m\ 



requires special treatment, 



pip is correct only for |m| 7^ |. The case 
because of the two degenerate roots in the upper component. For |m| = | 
the four linearly independent solutions behave asymptotically near the origin 
23 



like 



(Hp)\ 





1/P 



(22) 



We see that for any m there are always two solutions which are bounded near 
the origin and two others which are divergent. Since a smooth Hamiltonian 
can give rise to smooth eigenfunctions only |21| in the four dimensional space 
of solutions to the differential equation, there is a two dimensional subspace 
of admissible solutions, the ones which are well behaved at the origin. 



4.3 Solution to the ODE 

In this section we show that the solutions of Eq. (|]) that are regular at 
the origin, can be explicitly constructed in terms of certain Hypergeometric 
functions. 

Theorem 4.1 The solutions of ^ which are bounded at the origin are 
spanned by: 



(1)/ \ \ / n"^-k TP ( ■ ^ l,m 5,m.p''^ 



^ ^ \ 6+4m0-' 3li 372^ 3'6^ 3' 6*'' 

irP'^(n\ \ ( P™^^ Pr-S 3,m 7,m.p6' 
\ -[ VI^.^^A^ 3' 2 + T' 6 



•^m yP) I (2). X I 1 2 1 , m 7 , m.pS ' 



where 0-^3 (; a, b, c; x) are the generalized hypergeometric functions of the kind 

0-^3- 

Proof: Under scaling, p Xp, Eq. (H) transforms to 

-T- + + X m ^Mp) = 0. (24) 



dp2 pdp p2 \^ X-ip _ 



In particular, the equation is invariant under scaling by A, a cube root of 
unity, A = e^'^*/^. (Note that this feature is lost when one considers solutions 
of the equation for non-zero eigenvalue.) By an analog of Bloch theorem, 
the solution is a product of an eigenfunction of the scaling transformation 
and a periodic function under scaling by e^^^^^. p" is an eigenfunction of the 
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scaling transformation with eigenvalue A". Hence that Tm must be of the 
form p°'Q{p^). The indicial equation fixes a = m — 1/2. Qm{p^) is then an 
analytic function of its argument. 

The space of solutions regular of this kind is two dimensional and gives a 
representation of D^, the group of discrete rotations by 27r/3. Since the only 
complex irreducible representations |19| of are the complex numbers u, 
such that 00^ = 1, one can always find a basis J-'m\p) , J^m\p) such that 

J'^\e'-/'p)=u;,J^^\p). 

This condition fixes the solutions in Eq. (^) where Ui = e'^^^i"^-'^/'^)/^ and 

To relate Qm to hypergeometric functions we turn the two coupled second 
order equations into a scalar fourth order equation for each component. 
The equations obtained for the component and ip- can be written, with 
( = and D = in the form: 

The generalized hypergeometric function 0-^3 (; a, b, c; Q is defined by p^ : 

i7( h n ^ r(a)r(&)r(c) 

It is a matter of calculation to see that it satisfies the differential equation: 
{D(D + a - 1)(D + 6 - 1)(D + c - 1) - C} oFsi; a, b, c; () = 0. (27) 

Eq. (^) is a special case of this. Note, however, that we are not free to pick 
both ipY and (f^^ as the hypergeometric functions corresponding to Eq. (pT]). 
We can pick one, and then the other is determined by Eq. (^. 

To get jFrn"* we pick the upper component to be the hypergeometric func- 
tion that solves Eq. (|25|), i.e. 



3' 2 3 '6 3 ' 6^ 
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The lower component,*/?^^'', is determined by (H), which gives us the relation 



'^dp ydp J 



(28) 



With the identities 24, 26 



+ c- 1 ) oF3i;a,b,c;C) 



d 
dC 



oF3{;a,b,c; Q 



(c- l)oF3{;a,b,c- 1;C) 
1 



abc 



oF3(;a+l,6+l,c+l;C) 



we obtain 



, 1 ^, 5 ^ / 4 3 m 5 m p« 



(29) 



The second solution, J-'m\ is obtained by picking ip^^^ to be the hyperge- 



ometric solution to Eq. 



Namely, 



2 1 m 7 m p' 



+ 



+ 



3' 2 3 '6 3' 6 



(30) 



and with a relation 

.(2) 



,,,, , p'^+a ( d f d ^ 



ip^^\m;p) 
p^+l 



similar to ( pSj) one computes the upper component of the second solution: 



1 



12 + 8m 



5 3 m 7 m p^ 
3' 2 ^ 3"' 6 ^ 3"' 64 



4.4 The well behaved solutions 

We have seen that of the four dimensional family of solution of Eq.(H) there is 
a distinguished two dimensional family that is well behaved near the origin. 
We shall now show that there is a three dimensional family that is well 
behaved at infinity: 
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Theorem 4.2 • In the four dimensional space of solutions of Eq. 

there is a three dimensional family of solutions that vanish at infinity, 
and one dimensional subspace that diverges exponentially at infinity. 

• The solutions of (Qj for p >> 1 are (asymptotically) spanned by the 
four dimensional family : 

P-'I' exp 1 ^ . ^-3/4 0^3/2^ ^ 1 

P'"" %^'") ( -1 ) ' P'"" ( -1 ) • (3^) 

• The exponential blow up of Tm^^ of Eq. ^ is given by 

• The solution to Eq ^ that vanishes at the origin and at infinity is a 
multiple of 

Proof: H(,{p, ^,m) for Eq. (|]) is, for p » 1, 

HM.M^'^.{i S) (32, 

With eigenvalues ±p. It follows that the solution for large p reduces to the 
study of two uncoupeld equations: 

1 d \ 

±p)tp = 0. (33) 



dp"^ p dp 
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Since p is large, these can be solved by WKB to give the first part of the 
theorem. The blow up of the solutions at infinity can be obtained from the 
relation (omitting the exponentially decaying part) 

oF3(; a, 6, c; x) ^ ^-MF^F^^i {^^^'^ + 2 cos(4a;^ + 27r7)) , (34) 

with 7 = — °+^+^~^/^ _ This relation can be obtained by studying the asymp- 
totic behavior of the coefficients in the series of 0-^3- Alternatively, in [^, the 
asymptotic form of the generalized hypergeometric functions pFq is derived, 
and the formula given there reduces to ( P^ after substituting g = 0,p = 3, 
computing the summations and omitting the exponentially decaying part. 
iFiom. this the rest follows, as well as the proof of theorem |3.1| . ■ 
It remains to explain how eigenvectors are related to these well behaved 
solutions. The point is that the canonical differential equation approximates 
the eigenvalue equation only for r << 1, or, equivalently, for p « p~^^^. 
Consider an eigenfunction. Far from the crossing this eigenfunction can be 
approximated by a WKB solution, and it is clear that this WKB solution can 
be approximated by WKB solution of the canonical problem in the interval 
>> p >> 1. The component that blows up must have an exponentially 
small amplitude, of order exp(— g^), and, to leading order can be neglected 
near the crossing. 



5 Anomalous Mixing 

The basic and fundamental observation of the Born Oppenheimer theory is 
the emergence of the energy scale ^JJi in molecular spectra associated with 
vibrations. Perhaps the most interesting observation that results from the 
analysis of the Born Oppenheimer theory near crossing, is that emergence 
of the scale, p}^^^ associated with mixing at crossing, p^^^ is normally not a 
small number: In molecules p^^^ ~ .2. 

Classically, a uniform density on the energy shell, 5{p'^ + V [x) — E) implies 
that, in two dimensions, the spatial density is also uniform on the classically 
allowed region. The region where there is substantial mixing between the two 
electronic energy surfaces has linear dimensions that scale like p^/'^. For a 
crossing point in two dimension the volume characterizing the mixing there- 
fore scales like p^/"^. The semiclassical expectation is therefore that mixing 
near crossing should scale like the area /i^/^. 
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For isotropic crossing we found that the wave function has anomalously 
large amplitude in the mixing region for values of azymuthal quantum num- 
bers that are small compares to For these, from theorem |3.2| , the 
amplitude in the near zone scales like This implies that the total 
mixing weight scales like fi^^^. 

It would be interesting to have a more complete picture of the mixing 
near non-isotropic crossing and also for chaotic systems. 
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